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Executive  Summary 

This  preliminary  investigation  addresses  key  program  elements  for  sonar  sensing  in  a 
shallow'water  environment  before  development  of  a  full  Navy  proposal.  Such  an  investigation  is 
fundamental  to  the  establishment  of  bounds  on  possible  solutions  and  to  the  reduction  of  program 
uncertainty.  The  modeling  and  experimental  program  addressed  here  is  limited  to  two  issues — ^the 
potential  degradation  of  sonar  data  due  to  signal  masking  by  shallow-water  reverberation  and 
signal  loss  caused  by  extreme  platform  motions.  These  two  issues  were  selected  as  a  research 
focus  because  they  are  unique  to  the  shallow-water  environment  and  no  literature  exists  to  guide 
project  development. 

The  research  program  combined  theoretical  modeling  of  the  two  phenomena, 
experimental  validation  in  a  shallow-water  environment,  and  further  development  of  a  computer 
model  to  explore  parametric  sensitivity.  Models  simulating  the  signal  motion  loss  caused  by 
extreme  platform  motions  were  developed.  A  reverberation  model  for  estimating  the  average 
reverberation  intensity  in  the  shallow  water  environment  is  presented.  Computer  simulation 
results  show  that  the  performance  of  the  models  is  good.  These  results  are  compared  with  real 
sonar  data  from  an  initial  dock-side  test,  which  shows  good  agreement  with  the  theoretical 
prediction. 

As  a  result  of  the  shallow- water  experiments  and  acoustic  modeling,  we  conclude  that*  (1) 
Signal  motion  loss  is  not  the  dominant  factor  for  sonars  in  the  frequency  range  of  interest  (>200 
kHz);  rather,  a  high-quality  (velocity-aided)  inertial  navigation  and  attitude  system  will  be 
sufficient  to  correct  for  geometric  distortions  caused  by  platform  motion.  (2)  Although  surface 
reverberation  and  multipatii  noise  can  be  a  factor,  particularly  in  shadow-mode  imaging, 
reverberation  levels  are  rapidly  attenuated  at  the  frequencies  of  interest  and  beam  patterns  can  be 
manipulated  to  reject  most  interferences;  echo-mode  imaging  is  still  dominated  by  the  contrast 
between  target  strength  and  bottom  reverberation. 
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1.0  INTRODUCTION 

It  is  anticipated  that  future  wars  will  involve  third-world  scenarios  characterized  by 
special  operations/low-intensity  conflicts  (SOLIC)  with  amphibious  assaults  such  as  that  used  in 
the  Persian  Gulf  war.  Doctrine  calls  for  both  ingress  and  egress  assault  lanes,  each  approximately 
150-m  wide  and  up  to  2-km  long  for  a  regimental-sized  assault  force.  Current  above- water 
detonation  technology,  which  seeks  to  explode  mines  with  significant  overpressure,  has  been 
shown  to  be  mininudly  effective.  As  a  result,  it  is  expected  that  there  will  be  up  to  SO  mines  in 
areas  of  depth  less  than  10  ft,  which  will  be  cleared  using  such  direct-contact  techniques  as  towed 
sleds.  There  will  also  be  tqrproximately  10  to  IS  mines  in  the  shallow-water  area  between  10  and 
80  ft  of  depth.  These  mines  will  consist  of  both  older  metal-based,  contact  technology  and  newer 
composite-based  influence  mines  ranging  in  size  from  3-ft  diameter  in  the  deeper  water  to  6-in 
diameter  near  the  surf  zone.  The  mines  can  be  triggered  from  sensed  magnetic  or  electromagnetic 
fields,  sound,  overpressure,  or  contact. 

Special-warfare  personnel  are  currently  deployed  to  emplace  detonators  on  shallow-water 
mines  for  subsequent  destruction.  The  deployment  and  recovery  of  these  swimmers  is  extremely 
dangerous  because  of  premature  detonations,  malfunctioning  equipment,  and  enemy  defenders. 
MITRE’s  Autonomous  Systems  Laboratory  (ASL)  and  the  Deep  Submergence  Laboratory  (DSL) 
of  the  Woods  Hole  Oceanographic  Institution  (WHOI)  have  been  working  with  the  Navy  to 
develop  a  proposal  for  the  demonstration  of  small  underwater  vehicles,  supervised  by  humans 
operators  remote  from  the  hazardous  operation,  to  clear  mines  in  the  shallow-water  area.  We 
believe  that  such  semi-autonomous  underwater  vehicles  (SAUVs)  can  be  an  effective  asset  for 
shallow-water  mine  countermeasures  when  designed  and  equipped  with  appropriate  sensors  and 
actuators,  all  under  the  high-level  management  of  recently  developed  AI  techniques  for  situated 
reasoning. 

2.0  APPROACH 

In  considering  a  long-term  approach  to  the  problems  of  mine  countermeasures  in  a 
shallow- water  environment  by  an  SAUV,  we  generally  partition  our  approach  into  three 
components:  remote  sensing  and  machine  perception;  platform  engineering  and  closed-loop 
control;  and  reactive  tactical  control  and  mission  management  For  each  component  we  have 
focused  two  key  experiments  as  a  means  of  addressing  major  programmatic  uncertainties.  The 
results  of  these  experiments  will  be  used  as  initial  parameters  for  a  system-level  simulation  of  the 
larger  problem. 

A  full  treatment  of  the  problem  is  not  addressed  here.  Instead,  we  focus  on  sensing  issues 
germane  to  any  solution  for  detection,  classification,  and  task-relative  servoing,  whether 
accomplished  manually  or  semi-autonomously.  Successful  completion  of  such  mission 
components  will  be  largely  determined  by  the  quality  of  available  sonar  data.  The  two 
environmental  effects  expected  to  most  significantly  degrade  such  data  are  signal  masking  by 
shallow- water  reverberation  and  signal  loss  caused  by  extreme  platform  motions.  To  address 
these  two  issues  our  research  program  combined  theoretical  modeling,  experimental  validation  in 
a  shallow- water  environment,  and  further  development  of  a  computer  model  to  explore 
parametric  sensitivity. 
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Figure  1  Geometry  of  shallow- water  experiments  off  WHOI  dock,  showing  a  300-kHz 
forward-scan  sonar  mounted  on  a  remotely-operated  vehicle. 


2.1  Shallow- Water  Experiments 

To  support  theoretical  analysis  and  numerical  modeling,  experimental  testing  was 
conducted  in  a  shallow-water  environment  off  the  WHOI  dock  using  a  remotely  operated  vehicle 
(ROV).  Figure  1  depicts  the  (down-looking)  geometry  for  the  expoiments,  which  used  a  SOO-kHz 
forward-scan  sonar,  mechanically  rotated  with  a  precision  stepping  motor  (Fig.  2).  Although  the 
ROV  was  also  equipped  with  a  200-kHz  sidescan  sonar,  tight  scheduling  constraints  limited  die 
amount  of  time  available  for  testing,  and  efforts  focused  on  the  3(X)-kHz  system.  This  was  not  a 
severe  drawback  because  the  sonar’s  narrow  horizontal  beams  are  identical  to  a  typical  sidescan 
system  and  could  be  used  to  validate  both  side-  and  forward-look  geometries. 

The  forward-scan  sonar  is  a  member  of  DSL’s  family  of  phase-difference  systems,  which 
includes  120-kHz,  200-kHz,  and  3(X)-kHz  sonars  designed  for  seafloor  imaging  and  swath- 
bathymetric  mqiping.  The  sonars  were  developed  by  DSL  in  conjunction  with  the  Applied 
Physics  Laboratory  (APL)  of  the  University  of  Washington  and  Acoustic  Marine  Systems,  Inc.  of 
Redmond,  Washington. 
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Figure  2  300-kHz  scanning-sonar  geometry. 


The  system  has  been  fully  characterized  against  reference  hydrophones  at  the  APL 
acoustic  calibration  facility  to  quantify  source  level,  beam  pattern,  receiver  sensitivity,  and  system 
gain;  both  near-  and  far-field  responses  were  considered.  The  dual-receiver  design  also  provides 
phase  information  that  can  be  used  to  generate  high-resolution  swath  bathymetry;  system 
calibration  includes  phase-center  separation  and  acoustic-axis  offset 

Transmit  state,  repetition  rate,  pulse  length,  and  gain  of  the  systems  are  operator 
selectable.  Other  significant  features  include  digitizadon  at  the  receiver,  multiple  channels  of 
quadrature-detected  data,  and  a  high-bandwidth  optical-fiber  link  to  the  surface.  For  each  receiver 
channel  (two  per  beam),  a  surface-selectable  step  gain  is  applied  and  the  antialiased  signal  is 
quadrature  sampled  at  10  kHz.  The  digital  signal  is  Manchester  encoded  at  a  2-MHz  clock  rate 
and  transmitted  to  the  surface  over  a  fiber-optic  serial  link  (eight  16-bit  numbers  per  sampling 
interval).  At  the  surface,  the  serial  stream  is  decoded  and  digital  data  are  transferred  to  the  surface 
processor  and  merged  with  navigation  (pressure  depth,  acoustic  x-y-z)  and  attitude  (flux-gate 
compass,  gyro,  pitch,  roll)  data,  then  time  stamped  and  recorded  in  digital  form  on  8-mm  helical- 
scan  tape. 

3.0  SIGNAL  MOTION  LOSS  (SML) 

3.1  Background 

A  sidescan  sonar  signal  is  generally  displayed  and  evaluated  as  a  backscattered-intensity 
image.  For  such  images,  there  are  two  primary  types  of  distortion,  geometric  and  radiometric 
(intensity).  The  former  corresponds  to  a  discrepancy  between  the  relative  location  of  the  feature 
on  the  image  and  its  true  location  on  the  seafloor,  the  latter  is  a  deviation  from  the  ideal  linear 
relation  between  image  intensity  and  backscattering  strength  of  the  materials  on  the  seafloor. 
Various  kinds  of  distortions  can  be  identified  in  each  of  these  two  categories  [1],  [2]. 

Intensity  distortions  in  sidescan  sonar  images  arise  from  a  number  of  different  sources. 
Spherical  spreading  and  attenuation  in  water  cause  the  acoustical  intensity  to  decrease  as  sound 
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Hgure  3  Geometry  of  transmission  and  reception  for  a  backscattered  acoustic  signal. 
The  more  general  bistatic  case  encompasses  monostatic  operation. 


propagates  away  from  the  vehicle.  Another  factor  that  affects  the  signal  power  is  the  decrease  in 
backscattering  intensity  due  to  diminishing  grazing  angles  at  increasingly  longer  range.  Extensive 
research  has  been  carried  out  to  mitigate  the  problems  associated  with  such  distortions  [1]‘[6]. 
The  general  approach  is  to  apply  a  time-varying  gain  (TVG)  to  the  returned  signal,  allowing  the 
user  to  adjust  the  image  intensity  to  make  it  as  uniform  as  possible. 

However,  another  intensity  distortion,  which  corresponds  to  the  platform  motion  of  a 
vehicle,  remains  almost  untouched  by  previous  woiic.  The  platform  motion  could  cause 
unexpected  and  unpredictable  misalignment  of  the  transmitting  and  receiving  beam  patterns, 
resulting  in  a  drop  in  the  power  level  of  the  returned  signal.  Previous  investigations  have  not 
addressed  this  problem  for  three  reasons:  (a)  the  lack  of  precise  navigation  information  makes  it 
hard  to  evaluate  the  significance  of  the  problem,  (b)  the  traditional  use  of  sidescan  sonar  in  calmer 
deep  water  makes  it  less  significant  than  in  a  shallow-water  environment,  and  (c)  the  broader 
azimuth  beam  width  (at  typically  lower  frequencies)  makes  the  received  signal  level  relatively 
insensitive  to  platform  motion. 

3.2  Theoretical  Model 

Figure  3  [1]  illustrates  the  geometry  of  motion-loss  distortion.  At  time  to,  the  transmitter, 
at  position  x,,  emits  a  pulse  of  acoustic  energy  in  the  direction  a(6,  ^),  where  6  is  the  votical 
angle  and  is  the  azimuth  angle.  A  receiver  with  translational  and  rotational  components  of 
motion  receives  the  backscattered  acoustic  energy.  This  process  is  continuous,  starting  with  the 
transmission  of  the  sonar  pulse  until  it  reaches  its  maximum  detectable  range.  At  a  time  tj,  the 
sound  scattered  by  a  surface  patch  Xg  reaches  the  transducer  located  at  the  position  Xp  which 
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points  in  the  direction  Oj.  To  describe  the  motion  loss,  a  signal-motion-loss  (SML)  coefficknt  can 
be  defined  as: 


SML  it)  = 


1} 


B,  (0.  B,  (0. 0.  t)  BSM  (0,  t)  ds 


(0,^)BSW(0,^,/)dj 


(1) 


where  B^,  and  B((0, 4)  are  the  transmitting  and  receiving  beam  patterns,  and  is  the 
scattering  patch  on  the  range-surface  ellipse.  The  backscattering-strength  map  (BSM)  specifies 
the  backscattering  intensity  on  the  patch. 

In  general,  the  platform  has  six  degrees  of  freedom,  three  translational  and  three 
rotational.  However,  we  can  show  that  the  SML  caused  by  translational  motion  is  small  enough  to 
be  neglected.  Figure  4  compares  the  misalignment  of  the  beam  patterns  due  to  translational 
motion  with  that  due  to  rotational  motion.  Figure  4a  shows  the  effects  of  a  translational  motion  at 
10  m/s,  which  is  much  higher  than  normal  vehicle  operating  speeds.  However,  the  misalignment 
in  beam  patterns  can  hardly  be  seen  due  to  the  considerably  higher  sound  speed  (1500  nVs,  about 
two  orders  of  magnitude  greater  than  the  vehicle  speed). 

The  scenario  is  quite  different  for  the  rotational  case,  as  shown  in  Fig.  4b,  which 
corresponds  to  a  rotational  speed  10  deg/s.  The  percentage  of  the  misaligned  part  increases  with 
time,  assuming  a  constant  rotational  speed.  Although  the  model  developed  here  includes  the 
translational  motion  case,  calculations  show  that  the  SML  due  to  the  cross-track  translational 
motion  loss  is  less  than  1%  and  the  along-track  translational  motion  at  typical  vehicle  speeds 
causes  only  a  5%  drop  in  signal  levels.  In  the  following  discussion  we  therefore  neglect  the 
motion  loss  due  to  the  translational  motion  and  consider  only  the  rotational  case. 

3.2.1  One-Dimensional  Model 


We  begin  with  the  simplest  case:  a  one-dimensional  deterministic  motion  in  a 
homogenous  medium.  The  theoretical  beam  pattern  of  a  line  array  is 


B(0)  = 


rsin( 


nL 

X 


sin  (0) ) 


nL 

X 


sin  (0) 


2 


(2) 


where  L  is  the  array  length  and  X  is  acoustic  wavelength.  Figures  Sa  and  Sc  show  the  beam 
patterns  in  polar  coordinates  with  rotational  speeds  of  10  deg/s  and  20  deg/s  respectively;  Figs.  Sb 
and  Sd  show  their  counterparts  in  Cartesian  coordinates. 
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Figure  4  Shifted  beam  patterns  caused  by  translational  and  angular  platform  motions. 
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A  change  in  the  angle  causes  the  receiver  beam  pattern  to  move  along  the  axis.  The  SML 
is  determined  from  the  integral  of  the  overiapped  area.  The  one-dimensional  form  of  (1)  is 
derived  by  omitting  a  proportionality  constant  [7],  giving 


SML{u)  = 


where  u  is  the  angular  shift  due  to  platform  motion.  The  uniform  medium  assumption  implies 
BSM  is  constant  and  for  our  purposes  can  be  assumed  to  be  unity  without  loss  of  generality.  From 
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Figure  6  One-dimensional  signal  motion  loss  (SML)  for  constant  angular  velocity. 


Fig.  5  and  (3)  we  see  that  this  is  a  convolution  process.  With  simplifying  assumptions  we  derive 
the  analytical  expression  of  (3), 

SML  («)  =  — \-=  ( ^  [cos  ( {k%  -  ku)  -  cosku)  ]  ^  +  2kn 
7C  (2m  -  7t)  ' 

(4) 

-^sin  (Iku)  +^icos2ku  +  2)lni^^^J  ) 

Figure  6  shows  the  results  of  SML  computed  from  (4)  as  functions  of  time  and  rotational  speed 
using  a  constant  speed  for  each  curve,  illustrating  a  descending  trend  of  SML  with  increasing 
rotational  speed. 

3.22  Three-Dimensional  Model 

For  a  three-dimensional  problem,  difficulties  in  analytically  taking  the  double  integral  of 
(1)  lead  us  to  adopt  a  numerical  approach.  Figure  7  shows  the  3-D  representations  of  theoretical 
beam  patterns  assuming  a  reciangular  array  [4], 

JcL_  2 

£(6,  (|))  =  {  sine  (-^sinecos(|))  sine  (-^sinesin(j>)  }  (5) 
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Figure  7  Theoretical  transmit  and  receive  beam  patterns  for  a  rectangular  array. 


where  and  Ly  are  linear  dimensions  of  the  array  and  sinc(x)  =  sin(xyx.  The  rotated  receiver 
beam  pattern  Bf  can  be  calculated  by  rotating  the  beam  axis  using  three  rotation  matrices.  Figures 
7a  and  7b  are  the  transmitting  and  receiving  beam  patterns  as  in^cated;  Fig.  7c  shows  the  two 
beam  patterns  on  the  ‘^^me  graph.  Figure  7d  shows  the  result  of  multiplying  the  two  beam  patterns 
of  Figs.  7a  and  7b. 

The  three  angles  can  be  measured  relative  to  the  jc,  y,  and  z  axes  (Eulerian)  or  be  expressed 
as  pitch,  roll,  and  yaw,  depending  upon  angular  sensor  input  At  each  instant  the  vehicle  has  a 
yaw  rotation  speed  and  a  pitch  rotation  speed  Wj^.  The  total  effect  causes  the  vehicle  to  rotate 
around  an  axis  in  the  cross-track  plane  at  speed  W,  where  W^=  ^z-  The  axis  can  be  found 

as  shown  in  Fig.  8.  We  rotate  the  3-D  receiving  beam  pattern  around  this  axis  to  get  B^,  t). 
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Figure  8  Composition  of  platform  angular  velocity. 


By  taking  advantage  of  the  narrower  along-track  beam  of  sidescan  sonar,  we  can  simplify 
the  problem  to  two  dimensions.  The  scattering  field  on  the  ensonified  patch  can  now  be 
approximated  as  a  constant,  allowing  us  to  factor  out  the  BSM  function  from  the  integral  over  ^ 
and  define  a  new  variable,  the  vertical-motion-loss  coefficient,  as 


VMic(»,e)  =  ^  *  ■  ' - 


(6) 


The  advantage  of  such  a  definition  over  (1)  is  that  (6)  gives  pure  motion  effects  on  the  scattering 
field. 


Figures  9-11  show  computed  VMLC  using  (6).  The  roll  effect  is  simply  a  shift  of  VMLC 
in  6  and  is  not  presented  here.  Figure  9  shows  that  the  signal  loss  due  to  pitch  is  minimum  along 
the  beam  axis;  Fig.  10  shows  that  signal  loss  due  to  yaw  is  maximum  along  beam  axis.  Figure  1 1 
shows  the  composition  of  the  two  motions  giving  rise  to  a  shifted  motion-loss  pattern. 

In  practice,  it  is  very  difficult  to  measure  the  angle  of  rotation  at  each  instant  in  time 
corre^n^g  to  the  sampUng  interval.  Instead  we  can  measure  the  rotational  speed  when  the 
sound  pulse  is  transmitted  then  ^proximate  the  motion  as  linear  with  a  constant  q>eed.  We  can 
make  ^s  assumption  because  of  the  short  pulse  length  of  a  high-frequency  sonar  and  a  relatively 
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Figure  9  Vertical-motion-loss  coefficient  (VMLC)  for  a  pitch  rate  of  10  deg/s. 
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Figure  10  Votical-motion-loss  coefficient  (VMLC)  for  a  yaw  rate  of  10  deg/s. 
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Figure  1 1  Vertical-motion-loss  coefficient  (VMLQ  for  a  composite  angular  motion  rates 
in  pitch  and  yaw,  each  at  10  deg/s. 


Figure  12  For  a  sample  interval,  7,  angular  motion  is  assumed  constant 


slow  change  in  platform  motion,  as  shown  in  Fig.  12.  If  only  angular  measurements  are  available, 
instead  of  using  angular  rate  measurements,  we  may  use  an  approximate  form  for  the  angular  rate, 

(7) 

where  7  is  the  elq>sed  time  between  the  measurements  of  angular  motion  and  k  is  the  sample 
number. 
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3.3  Experimental  Results 

An  experimental  test  was  conducted  during  shallow-water  work  off  the  WHOI  dock  using 
a  stepper  motor  mounted  on  the  vehicle  to  rotate  a  forward-scan  sonar  (Fig.  2).  Because  of  the 
narrow  horizontal  beamwidth  of  the  forward-scan  sonar,  the  test  is  equivalent  to  induced  yaw 
motion  of  a  sidescan  (or  forward-looking)  sonar.  The  test  was  conducted  in  relatively  calm  water 
at  a  depth  where  surface-wave-induced  motion  was  negligible.  In  this  conHguration  the  vehicle 
was  static,  and  the  stepping  motor  could  be  used  to  simulate  induced  yaw  motions  at  a  measurable 
and  variable  rate. 

There  are  two  program  sequences  that  control  the  motor  motion  and  the  time  between 
successive  sonar  pings.  The  Hrst  program  commands  the  motor  to  remain  at  a  fixed  position  while 
the  sonar  is  pinging  and  collecting  data;  the  second  program  allows  the  CPU  to  synchronize  the 
pulse  transmission  with  sonar  scans  so  the  receive  beam  is  in  constant  motion  during  reception  of 
a  backscattered  pulse. 

The  first  data  set  corresponds  to  the  sonar  signal  without  motion  loss;  the  second  can  be 
thought  of  as  sonar  data  with  a  yaw-motion  loss.  We  can  then  compute  the  SML  by  taking  the 
ratio  of  these  two  data  sets.  By  normalizing  the  backscattering  cross  section  with  respect  to  area  (a 
range-dependent  function  of  beam  pattern  and  grazing  angle),  the  continuous  bottom  return  can 
be  considered  as  resulting  from  a  constant  target  strength  over  all  ranges. 

In  our  calculations,  we  simplify  the  problem  by  assuming  that  scattering  is  dominated  by 
direct  bottom  and  sea-surface  backscatter,  ignoring  scattering  contributions  from  other  ray  paths. 
For  bottom  backscattering,  we  use  Lambert's  law, 

5^  =  lOlog(p) +2Olog(sin0)  (8) 

where  5^  is  the  bottom  backscattering  strength,  p  is  the  bottom  scattering  constant,  and  6  is  the 
grazing  angle.  For  backscattering  from  the  sea  surface,  we  adopt  Marsh’s  empirical  model  [3], 

5,  =  -36-»-4Olog(tan0)  (9) 

where  Sj  stands  for  surface  scattering.  Figure  13  shows  a  typical  BSM  function  with  bottom  and 
surface  scattering.  Since  scattering  from  the  surface  (peak  around  vertical  angle  -SO  deg)  is  much 
smaller  than  that  from  the  bottom  (peak  around  vertical  angl<  50  deg),  we  intentionally 
exaggerate  the  former  by  about  40  dB. 

Figure  14  gives  a  comparison  of  the  predicted  SML  with  experiment  results  for  a  range  of 
about  75  m  (0.1  s  on  the  time  axis).  Although  the  comparison  is  not  exact,  both  data  and  the 
simulated  SML  show  similar  time  dependence,  which  suggests  that  our  motion-loss  model  is 
valid.  The  noisy  SML  is  largely  the  result  of  restricting  the  scanning  area  to  a  small  sector  so  that 
bottom  inhomogeneities  over  larger  areas  would  not  detract  from  results. 
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NORMALIZED  SCATTERING  FUNCTION 


BSM  function 


Azimuth  angle  *20  -100  Vertical  angle 

Figure  13  IVpical  backscatteiing-strength  map  (BSM)  showing  surface  and  bottom 
returns. 
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Figure  14  Comparison  of  predicted  and  experimental  results  for  signal  motion  loss 
(SML)  of  a  300-kHz  sonar  using  a  backscattered  bottom  return. 
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4.0  REVERBERATION 

4.1  Brief  Description  of  the  Acoustic-Reverberation  Model 

Difficult  problems  encountered  in  the  study  of  sound  propagation  in  a  range-dependent 
shallow-water  environment  are  mode  coupling  for  low-frequency  applications  and  waveguide 
multipath  for  high-frequency  applications.  So  far,  there  is  no  exact  solution  even  for  an  ideal 
wedge  (homogeneous  media  wiA  smoothed  boundaries)  with  a  penetrable  bottom.  Rough 
interfaces  (sea  surface  and  seabed)  and  inhomogeneities  within  both  the  water  column  and  the 
sediment  furtlwr  complicate  theoredcal  predictions.  In  addition,  bubble  scattering  near  the  sea 
surface  strongly  depends  on  statistics  of  the  rough  interface  and  bubbles.  The  statistics  are 
functions  of  time  and  such  environmental  factors  as  wind  q)eed. 

Thus,  an  acoustic  reverberation  model  in  a  realistic  shallow- water  waveguide  involves  at 
least  three  separate  issues:  wave  propagation  in  a  penetrable  waveguide,  scattering  from  rough 
interfaces,  and  the  volume  scattering  from  bubbles  and  sediment  inhomogeneities.  Each  of  these 
is  a  difficult  topic  in  underwater  acoustics.  A  general  acoustic-reverberation  model  should  take 
into  account  not  only  the  three  separate  aspects  mentioned  above  but  also  the  more  complicated 
interactions  among  them. 

It  is  impractical  to  obtain  accurate  predictions  for  such  a  problem,  at  least  at  the  current 
stage.  However,  it  is  possible  to  use  either  a  simplified  model,  an  empirical  model,  or  some 
combination  of  the  two  to  approximately  describe  acoustic  behaviors  in  a  shaUow-water 
environment  The  model  should  encompass  the  essential  physical  processes  and  can  be  used  to  at 
least  qualitatively  address  some  of  the  basic  acoustic  issues. 

It  is  not  our  intention  to  develop  any  “new”  theoretical  models  but  to  adapt  or  modify 
some  existing  models  to  fit  our  purpose.  The  acoustic  model  used  in  this  report  incorporate  the 
following  assumptions: 

a)  homogeneous  water  column; 

b)  perfect  reflection  from  sea  surface  (reflection  coefficient  is  -1);  and 

c)  constant  slope  with  small  roughness  compared  to  the  footprint 
Other  features  of  the  current  model  include: 

a)  accounts  for  multiple  bounces  (reflections),  but  only  the  first  scattering  either  from  the  sea 
surface  or  from  the  ocean  bottom; 

b)  iq>plies  an  empirical  model  to  surface  scattering; 

c)  iq)plies  a  composite-roughness  model  to  bottom  scattering,  which  includes  empirical 
models  for  volume  scattering  and  bottom  attenuation;  and 

d)  assumes  the  separability  between  incident  and  scattering  directions  for  scattering,  i.e., 

Fsca^^v  K)  -  where  is  the  backscattering  function  [8]. 
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4.2  Theoretical  Model 
42.1  Propagation 

Since  the  frequencies  of  interest  in  this  study  are  relatively  high  (wavelength  is  much 
smaller  than  the  scale  of  the  shallow-water  waveguide),  a  few  eigenrays  are  adequate  to  describe 
the  propagation  of  an  acoustic  wave.  The  geometry  of  our  problem  is  depicted  in  Fig.  ISa  for 
bottom  scattering  and  in  Fig.  ISb  for  surface  scattering,  where  only  four  rays  are  shown  for  each 
case.  Figure  16  illustrates  the  geometry  of  relevant  parameters  required  by  the  model  including 
propagation  distance,  outgoing  and  incoming  angles,  as  well  as  incident  and  scattering  angles. 
Variables  in  Fig.  16  are: 

Sq  source  location 

Rq  receiver  location 

5,  /**'  order  source  images 

/?,  order  receiver  images 

y  wedge  angle 

angular  coordinate  of  source 

angular  coordinate  of  receiver 

radial  coordinate  of  source  (OSq  =  OSi  =  OS^) 

r^o  radial  coordinate  of  receiver  (ORq  =  OR  i) 

angular  coordinate  of  source  image 

({>r  angular  coordinate  of  receiver  image 

Qr  reflection  grazing  angle 

6]  incident  grazing  angle  of  bottom  scattering 

02  scattering  grazing  angle  of  bottom  scattering 

hj  z  coordinate  of  source  image  (52^) 

hr  z  coordinate  of  receiver  image  (RiB) 

Yo  outgoing  angle  (transmit) 

Yi  incoming  angle  (receiving) 

A  horizontal  offset  (AB) 

horizontal  distance  between  source  image  and  the  scattering  point  (AQ) 

T}  distance  between  the  source  image  and  the  scattering  point 
(total  transmitting  distance  =  QSj) 

r2  distance  between  the  receiver  image  and  the  scattering  point 

(total  receiving  distance  =  QRi) 
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Figure  IS  Geometry  of  bottom  scattering  ray  paths  in  a  shallow-water  waveguide  for 
down-slope  propagation.  For  each  case  only  four  rays  associated  with  first- 
order  scattering  are  shown. 


Sea  surface 


Figure  15  Geometry  of  surface  scattering  ray  paths  in  a  shallow-water  waveguide  for 
down-slope  propagation.  For  each  case  only  four  rays  associated  with  first- 
order  scattering  are  shown. 
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Figure  16  Geometry  of  a  transmit  ray  showing  relevant  parameters. 


Let  ^  =  ri  +  r2  =  c/,  where  c  is  sound  speed  and  /  is  dme,  then 

ct  =  Jhl  +  ll  +  Jh^+(/^-A)\ 

where* 

hr  = 

A  =  r^ocos<^^-r^cos4>,. 

Substituting  variables  from  (11)  into  (10)  and  solving  for  /j,  we  have 

,  _  -B±Jb^-4AC 

^  2A  ’ 


(10) 


(11) 


(12) 
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where 


A2 

A  =  l-^ 

C  =  hf  +  A^-M^ 
R^-hl  +  hl  +  A^ 


(13) 


M  = 


2R 


The  plus  and  minus  signs  in  (12)  correspond  to  down-slope  and  up-slope  propagation, 
respectively.  We  can  now  express  other  variables  in  terms  of  Ij  as 


r2=  J/.2+(/,-A)2 


(14) 


and 


0^  =  ej+2y 


(15) 


The  ensonified  area,  A5,  is 


(16) 


where  X  is  the  pulse  length  and  is  the  equivalent  horizontal  beamwidth.  For  a  rectangular 
transducer,  we  use  the  relation  given  in  [4], 


lOlog^^  =  lOlog 


2kL 


-t-9.2. 


(17) 


where  L  is  transducer  length. 
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By  combining  the  above  results,  the  average  reverberation  intensity  can  be  expressed  in 
terms  of  the  scattering  coefficient,  ****  geometric  terms  given  by  (14)-(17) 

as 


/  = 


/=’(ei,e2)e‘ 


m'^2 


(18) 


where  Bf  is  the  reflection  grazing  angle,  Rf,  is  the  reflection  coefficient  for  a  planar  interface,  and  a 
is  the  attenuation  constant  A  more  genei^  expression  including  multiple  bounces  between  the 
sea  surface  and  the  bottom  can  be  obtained  by  replacing  the  reflection  coefficient  term  in  (18) 
appropriately. 

The  average  reverberadon  intensity  for  surface  scattering  has  the  same  form  as  (18).  For 
convenience,  in  the  surface-scattering  configuration  we  use  the  sea  surface  as  the  horizontal  axis. 
Minor  modifications  are  needed  to  obtain  the  corresponding  parameters  in  the  new  coordinate 
system. 

4J2  Scattering 

Because  of  inadequate  studies  in  the  area  of  bistatic  scattering  from  a  rough  interface,  a 
pseudobistatic  approach  [8]  is  used  in  our  model.  The  scattering  function,  Fsca/ki,  is  assumed 
to  have  an  approximate  form. 


PscaMi^K)  (19) 

where  Fi,j  is  the  backscattering  function;  A:,-  and  represent  the  incident  and  receiving  directions, 
respectively.  It  is  apparent  that  for  backscattering  (19)  gives  the  correct  answer.  Although  not 
accurate  to  describe  bistatic  scattering,  it  is  satisfactory  for  our  current  study,  as  in  most  cases  the 
eigenrays  of  direct  bottom  backscattering  contribute  most  to  the  total  scattering  field. 

For  scattering  from  die  sea  surface,  considering  that  transducers  are  mounted  with  a 
downward  tilt  in  most  cases,  the  beam  pattern  and  inherent  angular  dependence  of  surface 
scattering  make  scattering  from  the  sea  surface  insignificant  when  compared  with  the  total 
scattered  field,  except  at  the  instant  when  the  transmitted  signal  just  reaches  the  air-water 
interface,  which  results  in  a  strong  coherent-scattering  component  Furthermore,  the  analytical 
expressions  used  to  describe  scattering  firom  the  sea  surface  are  con^licated  because  of  the 
interactions  between  surface  and  volume  scattering.  Keeping  these  facts  in  mind,  we  use  a  single 
empirical  scattering  model  to  approximate  the  scattering  process  by  the  sea  surface.  The  nxxiel 
includes  wind  speed  and  bubble  scattering  effects.  In  our  computer  programs,  there  are  options 
for  several  different  empirical  models  [4]. 

For  bottom  scattering,  we  use  the  composite-roughness  model  to  incorporate  two-scale 
interface  roughness,  volume  scattering  of  the  sediment,  and  attenuation  in  the  water  and  sediment 
More  detailed  descriptions  of  this  model  can  be  found  in  [5],  [6]. 

Figures  17  and  18  show  the  time  series  of  reverberation  level  for  up-slope  and  down-slope 
propagations  with  a  source  of  unit  intensity.  In  these  figures  the  vertical  transducer  refuse  is 
shown  (i.e.,  for  a  normal  sidescan-sonar  geometry,  the  vehicle  would  be  moving  normal  to  the 
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REVERBERATION  LEVEL  (dB) 


Down  Slope  Propagation 


Figure  17  Time  series  of  reverberation  level  (RL)  for  down-slope  propagation  in  a 
shallow-water  waveguide.  The  spiky  curve  is  from  experimental  data;  the 
smoothed  curve  is  the  theoretical  prediction. 
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REVERBERATION  LEVEL  (dB) 


Figure  18  Tune  series  of  reverberation  level  (RL)  for  a  up-slope  propagation  in  a 

shallow-water  waveguide.  The  ^iky  curve  is  from  experimental  data;  the 
smoothed  curve  is  the  theoretical  pr^ction. 
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page).  The  0.2-s  tune  series  corresponds  to  a  range  of  ISO  m  (300-m  swath).  Hit  angle  denotes  the 
transducer  mounting  angle  from  the  horizontal  plane  (downward  is  positive).  The  spiky  curves 
are  experimental  data  collected  off  the  WHOI  dock;  the  smoodied  curves  are  corresponding 
theoretical  predictions. 

Some  parameters  used  in  computations  are  indicated  in  the  figures,  while  other  parameters 
related  to  bottom  and  subbottom  properties  and  required  in  the  simulations  are  consistent  with 
values  given  in  [6].  Since  the  system  was  not  fully  calibrated,  the  absolute  reverberation  level  was 
adjusted  to  obtain  best  fit  Although  the  absolute  level  is  a  free  parameter  in  Figs.  17  and  18,  the 
general  trends  of  the  reverberation  curves  agree  with  the  experimental  data  reasonably  well, 
which  further  confirms  our  acoustical  model. 

From  Figs.  17  and  18  we  can  see  that  as  expected,  the  reverberation  curve  for  up-slope 
propagation  decreases  more  slowly  than  for  down-slope  propagation.  With  the  same  tilt  angle  for 
both  up-slope  and  down-slope  propagation,  the  angle  between  the  incoming  wave  (resulting  from 
direct  bottom  backscatteiing)  and  the  beam  axis  is  smaller  for  up-slope  propagation  than  for 
down-slope  propagation.  A  smaller  angle  results  in  less  reduction  in  echo  level  due  to  the 
transducer  beam  pattern. 
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5.0  TARGET  DETECTABILITY 

In  this  section,  we  combine  two  models — motion  loss  and  acoustic  reverberation,  both 
discussed  separately  in  the  previous  sections — and  investigate  overall  sonar  performance  in  a 
shallow  water  environment  We  calculate  the  effective  veitical-beam-motion-loss  coefficient 
(VMLQ  defined  in  (6),  then  combine  it  with  the  2-D  vertical  receiving  beam  pattern  to  form  the 
effective  vertical  receiving  beam  pattern,  which  in  turn  is  used  by  the  reverberation-level  (RL) 
algorithm. 

We  use  an  echo-reverberation  ratio  as  a  reference  quantity  to  evaluate  the  sonar 
performance  in  terms  of  the  detectable  bottom  swath  width  (DBSW).  The  echo-reverberation 
ratio  represents  the  “signal”  (backscattering  cross  section  of  target)  to  “noise”  (spherical 
spreading,  attenuation,  and  beam-pattem-compensated  reverberation  intensity)  ratio,  which  is  a 
measure  of  target  detectability. 

We  use  the  DBSW  as  a  composite  measure  of  performance  for  the  echo-mode  detection  of 
both  midwater  and  bottom-proud  targets.  In  this  mode,  a  target  is  detectable  if  the  echo- 
reverberation  ratio  exceeds  a  given  detection  threshold  (DT).  This  measure  does  not  address 
shadow-mode  detection,  where  bottom  reverberation  levels  must  exceed  surface  reverberation 
plus  ambient-  and  system-noise  levels. 

Although  our  model  can  be  easily  used  to  predict  shadow-mode  performance,  our  efforts 
in  this  report  focus  on  echo-reverberation  ratio  for  the  purpose  of  demonstrating  modeling 
capabilities.  We  also  constrain  our  range  of  parameters  and  results  to  demonstrate  modeling 
C2q)abilities  rather  than  to  give  an  exhaustive  analysis  over  all  possible  sonar  systems  and 
environmental  conditions.  Our  aim  has  been  to  develop  a  tool  for  analysis  that  can  be  applied  to 
other  cases  in  support  of  a  more  focused  proposal  or  design  effort 

5.1  Echo-Reverberation  Ratio 

From  the  sonar  equation,  the  echo-reverberation  ratio  can  be  expressed  as  the  difference 
between  echo  level  (EL)  and  reverberation  level  [4], 

EL-RL  =  TS-  (S^-¥  lOlog A ) ,  (20) 

where  TS  represents  target  strength,  5^  denotes  scattering  strength,  and  A  is  the  equivalent 
scattering  area.  Both  terms  on  the  right  hand  side  include  influences  of  vehicle  motion  as 
discussed  in  Section  3.0.  Note  that  the  second  term  on  the  right-hand  side  (in  parenthesis) 
represents  a  pure  scattering  field. 

Figure  19  illustrates  such  a  pure  scattering  field  for  a  particular  case  where  a  random  noise 
background  has  been  added.  Since  the  noise  is  more  or  less  stationary  while  the  real  backscattered 
signal  from  the  surface  or  bottom  is  strongly  range  dependent,  at  larger  ranges  noise  will 
dominate  the  total  scattering  field.  Hence,  the  range  and  attenuation  corrections  in  computing  the 
scattering  strength  S^in  (20)  result  in  an  increased  scattering  field. 

Although  our  program  has  the  capability  to  model  targets  with  different  shapes  and 
material  properties  at  different  orientations,  for  simplicity  and  consistency  we  use  the  same  target 
throughout  our  analysis  since  the  echo-reverberation  ratio  is  a  function  of  target  strength  for  a 
given  geometry.  In  the  following  examples  we  used  a  prolate  spheroid  of  stainless  steel  with  a 
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Figure  19  2-D  image  of  echo  level  (EL)  in  a  shallow-water  waveguide.  Slope  is  0.3 

(wedge  angle  17.7  deg),  tilt  angle  is  10  deg,  3-dB  beamwidth  is  about  S3  deg. 
Source  and  receiver  are  located  at  (0,  -6.5). 
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broadside  incidence.  The  target  strength  was  calculated  using  a  nuxlified  high-pass  backscattering 
model  [9].  Although  the  high-pass  model  smooths  out  all  scattering  structures,  it  can  be 
considered  a  reasonable  model  as  long  as  an  averaging  process  is  used  [10].  More  detailed  studies 
of  target  scattering  involving  frequency  and  orientation  dependencies  can  be  fouiKl  elsewhere 
[11H14]. 

Figure  20  shows  the  image  plot  of  the  echo-reverberation  ratio,  the  contours  of  which  are 
shown  in  Fig.  21.  Figures  22-24  illustrate  echo-reverberation-ratio  contours  for  different  cases; 
the  open  circle  at  zero  range  indicates  source/teceiver  position  (backscatter).  Figure  22  shows  that 
a  narrower  beamwidth  gives  a  more  directional  detectable  area  (volume).  When  the  transducer 
looks  upward,  more  area  (volume)  in  the  upper  water  column  is  sensed  (Fig.  23);  looking 
downward,  more  area  (volume)  in  the  lower  water  column  is  sensed  (Hg.  24). 

To  study  the  dependence  of  the  target  detectability  more  efficiently,  we  define  DBSW  as 
the  swath  width  between  the  intersections  of  the  ocean  bottom  with  6-dB  echo-reverberation-ratio 
contour  lines  (horizontal  distance  between  a  and  b  in  Fig.  24).  We  ignore  the  region  near  nadir, 
which  results  from  a  null  of  the  transducer  beam  pattern.  Although  the  DBSW  is  defined  in  terms 
of  bottom  coverage,  this  is  only  a  convenient  (and  relevant)  measure  of  performance.  The  trends 
identified  in  following  results  apply  equally  well  to  midwater  targets. 

5.2  Results  and  Discussions 

All  results  shown  in  this  subsection  are  only  part  of  our  sensitivity  study  but  are 
representative  in  terms  of  DBSW.  Initially  there  were  11  candidate  variables  considered  important 
to  our  study;  transducer  depth;  tilt  angle  (transducer  mounting  angle);  transducer  length  and 
vndth,  which  control  beam  pattern;  bottom  slope;  bottom  depth;  acoustic  impedance  contrast 
across  the  seafloor,  which  affects  bottom  scattering  and  reflection;  wind  spe^,  which  affects 
surface  scattering;  and  vehicle  angular  velocities  (roll,  pitch,  and  yaw). 

Results  for  surface  scattering  are  not  shown  since,  according  to  early  modeling,  DBSW  is 
almost  unchanged  with  wind  speed  variation  within  a  reasonable  range,  and  DBSW  plots  are 
therefore  featureless.  Likewise,  DBSW  versus  vehicle  roll  motion  is  not  shown  because,  even 
though  it  affects  the  received  echo  level,  it  has  little  effect  on  DBSW.  This  is  because  both  target 
strength  and  reverberation  level  are  affected  approximately  proportionally  due  to  the  vehicle  roll 
motion.  Therefore,  the  following  analysis  focuses  on  the  remaining  nine  candidate  variables. 

Figures  25-40  are  representative  results  of  our  sensitivity  study  and  show  a  family  of  three 
DBSW  curves  plotted  against  one  of  the  nine  candidate  variables  (variable  1).  Each  individual 
curve  corresponds  to  a  chosen  value  of  one  of  the  eight  remaining  candidate  variables  (variable 
2).  The  other  seven  variables  are  held  constant  and  set  to  reasonable  default  values  listed  in  Table 
1.  In  the  foUowing  part  we  discuss  the  dependence  of  DBSW  on  all  nine  variables  in  more  detail 
and  explain  the  results  in  terms  of  the  underlying  physics. 
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Figure  20  2-D  image  of  echo-reverberation  ratio  in  a  shallow-water  waveguide.  All 

computational  parameters  are  the  same  as  in  Fig.  19.  The  target  is  a  stainless 
steel  prolate  spheroid  of  length  0.9  m  with  an  aspect  ratio  3:1.  Target  strength 
is  computed  on  the  assumption  of  broadside  incidence. 
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Figure  21  Echo-reverberation  ratio  contours  using  the  same  parameters  as  in  Fig.  20.  The 
open  circle  at  zero  range  indicates  the  source/receiver  location  (backscatter). 
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Figure  22  Echo-reverberation  ratio  contours  for  a  slope  of  0.3,  tilt  angle  of  10  deg,  and 
beamwidth  of  27  deg.  The  open  circle  at  zero  range  indicates  the  source/ 
receiver  location  (backscatter). 
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Figure  23  Echo-reverberation  ratio  contours  for  a  slope  of  0.1,  tilt  angle  of  -15  deg 
(looking  up).  The  open  circle  at  zero  range  indicates  the  source/receiver 
location  (backscatter). 
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Figure  24  Echo-reverberation  ratio  contours  for  a  slope  of  0.063  (wedge  angle  3.6  deg), 
tilt  angle  of  10  deg.  The  open  circle  at  zero  range  indicates  the  source/receiver 
location  (backscatter). 
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Table  1:  Default  Computation  Parameters 


Variable 

Value 

slope 

0.1 

total  depth* 

16.8  m 

transducer  depth 

6.S  m 

frequency 

300.0  kHz 

pulse  length 

200.0  ps 

transducer  length 

22.5  cm 

transducer  width 

5.0  mm 

tilt  angle^ 

10.0  deg 

sound  velocity  (water) 

1500.0  m/s 

sound  velocity  (sediment) 

1510.0  m/s 

water  density 

1.0  g/cm3 

sediment  density 

1.2  g/cir^ 

1.  depth  at  zero  range  (jcs=0) 

2.  transducer  mount  angle;  downward  is  positive 
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52.1  Transducer  Depth  (Figs.  25  and  36); 

In  general,  DBSW  increases  as  the  transducer  position  moves  closer  to  the  bottom.  This  feature 
can  be  understood  easily  by  realizing  that  as  the  transducer  approaches  the  bottom,  for  a  constant 
range  corresponding  to  a  particular  time,  the  backscattering  grazing  angle  decreases,  which 
produces  a  weaker  scattering  background.  Thus,  the  signal-to-noise  ratio  increases  and  results  in  a 
larger  DBSW. 

It  is  interesting  that  there  ate  discontinuities  observed  in  Fig.  36  (solid  line).  Such  a  discontinuity 
is  caused  by  the  sudden  separation  and  touching  between  the  6-dB  contour  lines  and  the  sea 
bottom.  This  phenomenon  can  be  seen  clearly  in  Figs.  21-23,  where  the  6-dB  contours  are  just 
above  the  ocean  bottom  on  the  down-slope  side.  There  are  many  places  where  such 
discontinuities  occur  in  other  figures. 

522  TiU  Angle  (Figs.  26, 27,  and  40): 

DBSW  increases  as  the  tilt  angle  increases.  It  is  obvious  that  more  “volume”  near  the  bottom  can 
be  ensonified  or  sensed  with  the  transducer  tilted  downwards. 

5.2.3  Transducer  Width — Vertical  Beam  Pattern  (Figs.  25, 26,  and  28-30): 

DBSW  increases  as  the  transducer  width  increases  (vertical  beamwidth  decreases).  Narrower 
vertical  beamwidth  results  in  less  scattering  from  the  ray  paths  going  out/coming  in  along  the 
directions  away  from  the  beam-axis.  Thus,  the  total  scattered  field  decreases  and  for  transducers 
tilted  facing  down  (positive  tilt  angle)  DBSW  increases. 

5.2.4  Transducer  Length — Horizontal  Beam  Pattern  (Figs.  28, 31,  and  37): 

Increase  in  transducer  length  is  equivalent  to  a  narrower  horizontal  beamwidth  and  a  smaller 
ensonified  area.  This  eventually  manifests  itself  as  less  scattering  from  both  the  bottom  and  the 
sea  surface,  resulting  in  a  larger  DBSW. 

5.2.5  Vehicle  Motion: 

There  are  three  different  results  with  respect  to  different  vehicle  motions: 

(a)  DBSW  decreases  as  yaw  speed  increases  (Figs.  29,  34,  and  35); 

(b)  DBSW  increases  as  pitch  speed  increases  (Figs.  31,  32,  and  38);  and 

(c)  DBSW  is  insensitive  to  roll  speed  (not  shown). 

For  case  (a),  with  the  given  transducer  beam  patterns,  the  larger  the  speed  of  motion  the  smaller 
the  overlapped  volume  between  transmitting  and  receiving  beam  patterns.  Since  at  any  particular 
time  the  motion  loss  is  maximum  in  the  direction  of  the  beam  axis  (Fig.  10),  the  reduction  in 
scattered  field  comes  mainly  from  the  direct  bottom  backscattering,  whose  path  is  around  the 
beam  axis.  At  the  same  time,  since  we  detect  swath  width  at  the  bottom,  the  signal  scattered  from 
the  taiget  (assumed  to  be  located  at  the  bottom)  decreases  at  the  same  rate  as  that  from  the  direct 
bottom  backscattering,  while  scattering  along  other  ray  paths  decreases  less  than  that  from  the  ray 
path  along  the  beam  axis  (i.e.,  due  to  vehicle  motion  the  percentage  reduction  of  the  total 
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scattoed  field  without  the  scatterer  is  less  than  the  percentage  reduction  in  the  target  strength). 
This  unbalanced  reduction  causes  a  lower  signal-to-noise  ratio  and  results  in  a  smaller  DBSW. 

Another  interesting  phenomenon  is  that  a  larger  target  gives  a  stronger  dependence  of  DBSW  on 
yaw  speed.  This  can  be  seen  by  comparing  Fig.  34  with  Fig.  35,  where  the  volume  of  the  target 
used  in  Fig.  35  is  about  16  times  as  large  as  that  used  in  Fig.  17.  When  a  larger  target  is  used,  the 
DBSW  becomes  larger  and  its  detectable  boundary  corresponds  to  a  later  tinw  and  a  lower 
reverberation  level.  Hence,  the  noise  contribution  to  the  total  scattered  field  becomes  more 
important  Because  the  noise  field  is  independent  of  motion,  the  reverberation  (noise)  background 
remains  unchanged  or  is  slightly  reduced.  On  the  other  hand,  the  backscattered  signal  from  the 
target  is  also  reduced  due  to  the  vehicle  motion,  but  as  a  percentage  it  is  more  than  that  for  the 
reverberation  background.  Hence,  the  signal-to-noise  ratio  decreases  faster  than  in  the  smaller 
target  case. 

For  case  (b)  (pitch),  the  motion  loss  is  minimum  along  beam  axis  (Fig.  9).  Following  the 
discussions  for  case  (a),  we  can  conclude  that  the  percentage  reduction  of  total  scattered  field  in 
the  waveguide  without  the  scatterer  is  larger  than  that  Ifiom  a  single  target.  Thus,  a  larger  DBSW 
is  expected. 

For  case  (c),  although  roll  changes  the  received  signal  level,  it  also  simultaneously  changes  the 
received  signal  from  the  target.  In  addition,  since  the  vertical  beamwidth  is  much  larger  than  the 
horizontal,  in  contrast  with  case  (4),  the  equivalent  beam  pattern  due  to  roll  doesn’t  change  very 
much  unless  the  speed  of  motion  is  very  large.  Therefore  DBSW  remains  virtually  unchanged. 

5.2.6  Bottom  Scattering  (Figs.  33  and  39): 

DBSW  in  general  decreases  as  acoustic-impedance  contrast  increases,  but  the  dependency  is 
weak  as  the  bottom  scattering  is  relatively  insensitive  to  the  impedance  contrast  at  the  seafloor 
interface  for  small  grazing  angles. 

5.2.7  Bottom  Slope  (Figs.  27,  30,  32,  and  34-39): 

DBSW  decreases  as  slope  increases  (Note  that  the  actual  slope  in  negative.  For  simplicity,  we  use 
absolute  value  of  slope  throughout  this  report).  As  the  slope  increases,  backscattering  grazing 
angle  increases  for  up-slope  propagation  and  decreases  for  down-slope  propagation.  An  increase 
or  a  decrease  in  grazing  angle  causes  a  decrease  or  an  increase,  respectively,  in  the  backscattering 
field.  As  discussed  in  Section  4.2.2,  reverberation  for  up- slope  propagation  decreases  more  slowly 
than  for  down-slope  propagation;  hence,  the  change  in  scattering  field  due  to  slope  on  the  up- 
slope  propagation  side  appears  to  be  more  important  than  that  on  the  down-slope  propagation 
side.  Thus,  it  is  expected  that  the  amount  of  decreased  DBSW  on  the  up-slope  propagation  side  is 
larger  than  the  amount  of  increased  DBSW  on  the  down-slope  propagation  side. 

52.8  Bottom  Depth  (Figs.  33  and  40): 

It  is  observed  that  DBSW  decreases  as  bottom  becomes  deeper.  As  bottom  depth  increases,  for  a 
fixed  range,  the  backscattering  grazing  angle  at  the  bottom  is  larger  than  that  with  a  shallower 
bottom.  The  increased  backscattering  field  due  to  a  larger  grazing  angle  results  in  a  smaller 
DBSW. 
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Figure  25  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  depth  and 
transducer  width. 
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Figure  26  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  tilt  angle 
and  transducer  width. 
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Figure  27 


Figure  28 


Detectable  bottom  swath  width  (DBS W)  as  a  function  of  transducer  tilt  an  ele 
and  bottom  slope.  ^ 


Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  width  and 
transducer  length. 
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Figure  29  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  width  and 
yaw  rate. 


Figure  30  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  width  and 
bottom  slope. 
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Figure  3 1  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  transducer  length  and 
pitch  rate. 
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Figure  32  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  pitch  rate  and  bottom 
slope. 
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Figure  33  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  sediment-water 
impedance  contrast  and  normalized  depth  (depth/Xdepth). 


Figure  34  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  yaw  rate  and  bottom 
slope. 
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YAW  Speed  (deg/s) 


Figure  35  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  yaw  rate  and  bottom 
slope  Gaser  target). 


Figure  36  Detectable  bottom  swath  width  G5BSW)  as  a  function  of  bottom  slope  and 
normalized  transducer  depth  (Xdepth/depth). 
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Figure  37  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  bottom  slope  and 
transducer  length. 


Figure  38  Detectable  bonom  swath  width  (DBSW)  as  a  function  of  bottom  slope  and 
pitch  rate. 
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Figure  39  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  bottom  slope  and 
sediment-water  impedance  contrast 


Figure  40  Detectable  bottom  swath  width  (DBSW)  as  a  function  of  depth  and  transducer 
tilt  angle. 
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5.3  Suinm£.ry  of  Sensitivity  Analysis 

To  conclude  our  sensitivity  studies  and  to  offer  a  more  intuitive  overview  of  the  above 
discussions  on  the  dependency  (sensitivity)  of  DBSW  on  the  different  variables.  Table  2 
summarizes  the  results.  Variables  ^m  left  to  right  (symmetric  to  those  from  top  to  bottom)  are 
arranged  according  to  their  relative  influences  on  DBSW.  In  other  words,  transducer  depth  causes 
DBSW  to  increase  the  most  while  bottom  depth  causes  DBSW  to  increase  the  least  (decrease  the 
most). 

The  value  of  subscript  of  I  (increase)  or  D  (decrease)  reflects  the  degree  of  the  change  in 
DBSW,  defined  as: 

0:  0<  5svkr<=  10(m) 

1: 10<  6s^<=20  (m) 

2:  5s^  >  20  (m) 

where  figw  represents  the  change  in  DBSW.  Mixed  behavior  (denoted  by  M)  means  there  was  no 
apparent  or  consistent  increasing  or  decreasing  trend.  In  each  rectangle  in  the  table,  there  are  two 
values;  the  one  above  the  diagonal  dashed  line  describes  the  dependence  of  DBSW  on  the 
variable  in  the  same  column  (variable  1)  while  the  one  below  the  diagonal  describes  the 
dependence  of  DBSW  on  the  variable  in  the  same  row  (variable  2).  We  note  here  that  the  decision 
of  whether  DBSW  is  an  increase  (denoted  by  I)  or  a  decrease  (denoted  by  D)  in  Table  2  is  based 
on  a  visual  average  over  the  three  curves  in  each  plot  It  is,  thus,  a  qualitative  measure  of  how 
DBSW  changes  in  terms  of  different  variables. 
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TABLE  2.  Performance  Sensitivity 


I:  Increase  in  DBSW 
D:  Decrease  in  DBSW 

M:  Mixed  behavior  (no 
apparent  increasing 
or  decreasing  trend) 
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6.0  CONCLUSIONS 

In  this  study,  we  have  investigated  the  overall  sonar  performance  in  terms  of  the 
detectability  of  a  tai^et  We  have  included  the  effects  of:  (a)  vehicle  motion  resulting  in  a  motion 
loss  in  most  cases;  and  (b)  reverberation  in  a  shallow  water  waveguide  with  small  roughness 
boundaries  and  a  penetrable  bottom.  Although  the  naodel  is  only  approximate,  the  results 
qualitatively  provide  us  with  parameters  for  defining  sonar  performance  under  different 
circumstances. 

In  general,  DBSW  has  stronger  dependence  on  the  geometric  configuration  of  the  sonar 
system,  such  as  transducer  position  and  tilt  angle,  than  on  such  a  physical  properties  of  the 
environment  as  acoustic  impedance.  In  another  words,  it  is  possible  to  optimize  the  system 
configuration  to  obtain  the  best  sonar  performance. 

It  is  surprising  that  although  vehicle  motion  can  greatly  affect  the  echo  level,  in  most  cases 
DBSW  is  relatively  insensitive  to  vehicle  motion  and  the  results  should  be  tested  whenever 
possible.  The  weak  dependence  of  DBSW  on  bottom  acoustic  impedance  contrast  is  also  quite 
surprising. 

It  is  worth  pointing  out  that  since  we  used  a  simple  quantity,  DBSW,  as  a  standard  in 
evaluating  sonar  performance  and  because  of  the  discontinuities  in  DBSW  observed  in  many 
plots,  DBSW  may  not  reflect  the  true  detectable  volume  of  interest  in  some  circumstances. 
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issues-the  potential  degradation  of  sonar  data  due  to  signal  masking  by  shallow-water  reverberation  and  signal  loss  caused  by 
extreme  platform  motions.  The  research  program  combines  theoretical  analysis,  experimental  validation  in  a  shallow-water 
environment,  and  development  of  a  computer  model  to  explore  parametric  sensitivity.  Results  from  an  initial  dock-side  test  show 
good  agreement  with  the  theoretical  predictions.  From  the  shallow- water  experiments  and  acoustic  modeling  we  conclude  that:  (1 ) 
Signal  motion  loss  can  influence  the  reverberation  level  significantly  but  is  not  the  dominant  factor  in  target  detection  for  sonars 
in  the  frequency  range  of  interest  (>200  kHz);  a  high-quality  (velocity-aided)  inertial  navigation  and  attitude  system  will  be 
sufficient  to  correct  for  geometric  distortions  caused  by  platform  motion.  (2)  Although  surface  reverberation  and  multipath  noise 
can  be  a  factor,  particularly  in  shadow-mode  imaging,  reverberation  levels  are  rapidly  attenuated  at  the  frequencies  of  interest  and 
beam  patterns  can  be  manipulated  to  reject  most  interferences;  echo-mode  imaging  is  still  dominated  by  the  contrast  between  target 
strength  and  bottom  reverberation. 
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